arXiv:1507.01785vl [quant-ph] 7 Jul 2015 


Dynamical moments reveal a topological quantum transition in a photonic quantum 

walk 


Filippo Cardano/ Maria Maffei/ Francesco Massa/ ’[^ Bruno Piccirillo/’^ Corrado de Lisio/’^ 
Giulio De Filippis/’^ Vittorio Cataudella/’^ Enrico Santamato/ and Lorenzo Marrucci^’^ 

^ Dipartimento di Fisica, Universitd di Napoli Federico II, 

Complesso Universitario di Monte Sant’Angelo, Napoli, Italy 
^ CNR-SPIN, Complesso Universitario di Monte Sant’Angelo, Napoli, Italy 


Many phenomena in solid-state physics can 
be nnderstood in terms of their topological 
properties mm- Recently, controlled protocols 
of quantum walks are proving to be effective 
simulators of such phenomena [SHB]. Here we 
report the realization of a photonic quantum 
walk showing both the trivial and the non-trivial 
topologies associated with chiral symmetry in 
one-dimensional periodic systems, as in the Su- 
Schrieffer-Heeger model of polyacetylene We 
find that the probability distribution moments 
of the walker position after many steps behave 
differently in the two topological phases and 
can be used as direct indicators of the quantum 
transition: while varying a control parameter, 
these moments exhibit a slope discontinuity at 
the transition point, and remain constant in the 
non-trivial phase. Extending this approach to 
higher dimensions, different topological classes, 
and other typologies of quantum phases may 
offer new general instruments for investigating 
quantum transitions in such complex systems. 

The presence of topological order in matter is respon¬ 
sible for fundamental phenomena, such as the fractional 
and integer quantum Hall effects 011 or the existence 
of protected surface states in topological insulators pT] , 
Non-trivial topological phases are related to symmetries 
and can be characterized by specific topological invari¬ 
ants, defined in terms of the energy band eigenstates. 
The properties of these topological phases can be con¬ 
veniently studied using quantum simulators 0IIIHI1. 
In this context, quantum walks (QW) are emerging as 
a versatile tool [4], since this simple quantum dynamics 
can realize all topological phases occurring in one- and 
two-dimensional systems of non interacting particles; 
as a first application of this concept, the formation of 
topologically protected bound states was observed in 
a photonic QW [S]. In this paper we report the study 
of a novel QW process that exhibits both the trivial 
and non-trivial topological phases associated with chiral 
symmetry in a unidimensional (ID) bipartite lattice. A 
remarkable example of a physical system having identi¬ 
cal symmetries and topological features is given by the 
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polyacetylene chain, whose mobile electron dynamics is 
well described by the Su-Schrieffer-Heeger model (SSH) 
0. SSH involves electron hopping between adjacent 
sites A and B, which may belong to the same cell or not, 
as shown in Fig. [It- Varying the relative strength of 
these two inter-site couplings, one predicts the existence 
of two distinct phases, marked by different values of 
a suitable topological invariant W. This topological 
classification is possible thanks to the presence of the 
chiral (or sublattice) symmetry associated with the 
interchange of the sublattices A and B mam] (see 
Supplementary Information [SI] for further details). By 
studying our QW system, we found that the probability 
distribution moments for the particle position after a 
long temporal evolution show a different asymptotic 
behavior in the two topological phases: varying an ex¬ 
ternal control parameter, these moments exhibit a slope 
discontinuity at the quantum transition, and remain 
locked to a constant value in the non-trivial topological 
phase. These effects occur also in the SSH model and 
hence seem to be fairly general, at least in the class of 
chiral-symmetric systems. We simulated experimentally 
such phenomena and propose a theoretical interpretation 
in terms of the dispersion relations and of the geometric 
features of the system eigenstates. Remarkably, our 
analysis takes into account bulk dynamics only, while not 
considering physical effects manifesting at the edges; the 
latter are the typical signature of topological systems. 
In the same spirit, topological phases have been re¬ 
cently detected in a non-Hermitian quantum walk 0112 ]. 

A ID quantum walk is the discrete evolution of a 
particle (the walker) on a ID lattice [T7]. At each 
step, the walker can move to the two nearest-neighbour 
sites; the direction of the shift is determined by the 
configuration of an internal 2-state quantum system (the 
coin). Between consecutive steps, the coin state under¬ 
goes a rotation, mimicking the coin toss characterizing 
the familiar classical random walk. QWs have been 
implemented in a variety of physical architectures [ 18 j . 
In our photonic platform dl, the walker is encoded 
in the orbital angular momentum (0AM) of light |2()| : 
discrete positions on the lattice are associated with 
states jm), where m is an integer, describing a photon 
carrying mh of 0AM along its propagation axis. The 
coin is encoded in the polarization degree of freedom: 
vectors \L) and \R), representing left and right circular 
polarizations, respectively, are the two internal states 
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that determine opposite shift directions in the lattice. 
The quantum state of the photon after n steps of QW 
is denoted as IV'n) = where IV'o) is the initial 

(input) state and Uq is the unitary evolution operator 
of a single QW step (its formal expression is in Eq. 
Bll). This evolution operator is realized with two 


optical elements: a quarter- wave plate (QWP) and a 
q-plate (QP) (Eqs. |B11|B13 ). A g-plate is essentially 
a liquid-crystal birefringent cell having the optic axis 
arranged in a singular pattern [2T], with topological 
charge q (in our case q = 1/2). This patterned bire¬ 
fringence gives rise to an engineered optical spin-orbit 
coupling that induces the polarization-controlled shift 
of 0AM. As specified in Eq. |B13[ besides the charge q, 
the action of this device is determined by the value of 
the optical retardation 5, which can be tuned with an 
applied electric field [22]. A topological characterization 
of QWs with chiral symmetry may be defined in terms 
of the eigenstates of the operator C/q (see the SI) which, 
due to the lattice translation symmetry, are associated 
with two quasi-energy bands parametrized by a quasi - 
momentum k, as shown in Fig. (see also Eq. Al 


As for momentum k, a periodic quasi-energy variab. 
E replaces here the usual energy as a consequence of 
time (= step number) being discretized. The coin part 
of the band eigenstates is denoted as \(j)s,s{k)), where 
s £ {1,2} is the band index. Conveniently, these states 
can be represented as points on the Poincare sphere for 
light polarization, individuated by a three-dimensional 
unit vector that we refer to as ns{k) (see Eq. B15|. 
Chiral symmetry constrains vectors ns{k) to lie on a 
great circle of the sphere. When k varies throughout 
the Brillouin zone {—tt, tt}, the number of closed loops 
(winding number) described by ns{k) is the topological 
invariant W. As shown in Fig. [^, the value of 5 in the 
range {0, 27r} determines the existence of two different 
topological phases: a non-trivial phase with W = 2q = 1 
occurs when 5i < 5 < (52, while W = 0 (trivial phase) 
in the remaining region. In Fig. & it can be noted that 
when i5 = (5i and S = 62 the two bands are touching 
and the dispersion is locally linear at fc = 0 and fc = tt, 
respectively, as typically occurs for quantum transitions 
between topological phases | 2 ]. Interestingly, c ompa ring 
the expressio ns of group velocity Vs = dEs/dk { A21 and 
vector ns{k) (B15), we can observe that Uz = —Uy = Vs- 
Thus, in the non-trivial phase, the maximum of the 
group velocity is independent of 6 (see Fig. |^in the SI 
for further details). A recently introduced QW protocol, 
that is the Split-Step QW, has very similar properties 


The features of the energy bands and associated eigen¬ 
states have profound consequences on the system dynam¬ 
ics, which shows marked differences in the two topolog¬ 
ical phases [3lll5|. Here we characterize such different 
behavior through the analysis of the moments of the 
probability distribution P{m) associated with the walker 
position, defined as Mj = PiiTi). In particu¬ 


lar, we consider a photon starting its walk in the po¬ 
sition m = 0 with an arbitrary polarization, that is 
IV'o) = |0) ® |</o), where \(j)o) is a generic state in the 
coin Hilbert space. Simulations of the QW evolution in 
the large step-number limit show that the moments Mj 
assume a constant value (independent of d) in the non¬ 
trivial phase, and undergo abrupt variations at the phase 
transitions, that is at ^ = {i5i,52}. In the infinite-steps- 
limit, these moments have simple asymptotic expressions 
in terms of the energy band dispersion relations; in par¬ 
ticular, those relative to the first and second moments 
Ml and M 2 are the following (see the SI for a proof): 

Mi/n= {sy - Sz)L{5) + 0{l/n) (1) 

M 2 /n^ = L{ 6 ) + 0{l/n^) (2) 


where Si = ((/olCTilf/o), with i £ {x,y^z'\, are the expec¬ 
tation values of the Pauli operators Ui in the coin space, 
calculated with respect to coin initial state |(/o)- Inter¬ 
estingly, we observe that M 2 is independent of the initial 
coin state. The quantity L{ 6 ) appearing in Eqs. 1|2 


IS 


equal to the square of the group velocity Vs, and hence 
to the square of Uy, integrated in momentum space over 
the Brillouin zone: 


m = f 

j —• 
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For our QW process, th e in tegral reported in Eq. 
admits a closed form (Eq. C9). In Fig. [3^, the quantity 
L{5) is plotted as a function of (5; as a consequence of 
the discontinuity present in the group velocity, L{ 6 ) is a 
piecewise function, and manifests abrupt variations at (5i 
and 82 - In the non-trivial phase, it is equal to a constant; 
this result can be qualitatively understood by looking 
at how the group velocity dispersion is modified by a 
change of 5 (Fig. [^. On the other hand, in the trivial 
phase L{5) increases or decreases in the two regions 
{0,7r/2} and {37r/2, 27r}, respectively. A similar behavior 
occurs in the SSH model (see Fig. and hence seems 
to be a general feature of the non-trivial topology. In a 
QW with a finite number of steps, statistical moments 
Ml and M 2 have a continuous behavior, converging to 
that given by Eqs. m asymptotically as n —>■ 00 . For 
M 2 , this convergence is rapid and visible for values of n 
that are small enough to be achieved in an experimental 
simulation, whereas for Mi such process is much slower 
(see Fig. [^. Thus, as a figure of merit for the quantum 

, - //^2\ 

transition, we chose to analyze VM 2 /rV- = —. In 

particular, we implemented a 6-step QW and measured 
the latter quantity when varying the external parameter 
(5, so as to simulate the topological quantum transition. 
The layout of the apparatus is shown in Fig. A 
standard “heralded single photon source” is realized at 
the input of the QW system, as discussed for instance 
in Ref. [19]. After the initial state preparation (see the 
caption of Fig. for more details) the photon undergoes 
the QW evolution and, at the exit, is analyzed in both 
polarization and 0AM so as to determine the output 
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probabilities and the associated moments. Letting 6 vary 
in the range {0,7r} with steps of 7r/16, we determined 
the corresponding probability distribution of the walker 
after a 6-step QW; a typical example is reported in Fig. 
[^. In Fig. 1^ we plot the measured values for 
as a function of 6 . In good agreement with quantum 
predictions, the data are close to the asymptotic limit 
reported in Eq. the emergence of an abrupt slope 
variation at d = 7r/2 can be appreciated and it is evident 
that in the non-trivial phase M 2 is constant. The 
observed discontinuity is a signature of the underlying 
quantum transition, while the existence of a constant 
value for M 2 in the non-trivial phase reflects the features 
of the group velocity dispersion, and of the geometry 
and topology of the system eigenstates. Preparing the 
initial state of the coin (m=0 for the walker) in two 
non-orthogonal polarizations we also verified that M 2 
(unlike Mi) is independent of the coin initial conditions, 
as expected from Eq. [see Figs. c)]. This aspect was 
further investigated by measuring i/M^/n for several 
initial polarizations, corresponding to specific points 
along a meridian of the Poincare sphere, and repeating 
the experiment for two values of S in each topological 
sector. The final data, reported in Fig. [^, match well 
the predicted results and clearly show that the value 
of the second-order moment is independent of the coin 
initial state. 


In conclusion, we have shown that a signature of quan¬ 
tum transitions between distinct topological phases is 
present in the behavior of suitable dynamical observables 
that are easily accessed experimentally. In the context of 
ID systems with chiral symmetry, such as the SSH model, 
we propose that the statistical moments of the particle 
probability distribution in space in the large step-number 
limit provide a convenient choice of such observables. We 
experimentally validated this proposal by simulating this 
topological environment within a specific photonic QW 
architecture. The measured asymptotic moments as a 
function of a control parameter remain constant in the 
non-trivial topological phase and show slope discontinu¬ 
ities at the phase change. We also proved that these sta¬ 
tistical quantities can be simply related to properties of 
the energy bands and in particular are directly linked to 
the square of the group velocity, integrated over the Bril- 
louin zone. The latter, in turn, can be traced back to the 
geometric features of the system eigenstates (although 
not just to the associated topological invariant), thus 
providing insight on the general features leading to the 
emergence of such phenomena. In prospect, our approach 
based on dynamical moments could be applicable to the 
investigation of other classes of quantum transitions and 
to topological phases associated with more complex sym¬ 
metries or with a higher dimensionality, thus helping to 
shed new light on the physics of topological phenomena 
and on their simulation in suitable experimental archi¬ 
tectures. 
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FIG. 1. Topological characterization of the QW process and SSH model of the polyacetylene chain, a) Scheme of the poly¬ 
acetylene chain in the SSH model. The latter is described as a sequence of two species of carbon atoms, that we denote as 
A and B, which are linked by alternating single and double bonds. The unit cell is made of a pair of A and B atoms, whose 
repetition gives rise to the discrete lattice for the electron dynamics. The couplings between A and B sites in the same cell 
or in neighbor ones are represented by dashed and continues arrows, respectively. The ratio of the corresponding transition 
amplitudes determines the topological features of the system, b-c) Dispersion relation for the quasi-energy Es{k) and the group 
velocity Vs{k) of the Bloch bands in the QW. The dispersion curves for the two bands depend on the external parameter S; 
being the temporal coordinate a discrete variable, the quasi-energy is defined in a Brillouin zone {—tt, tt}. Here we report few 
examples; colors of the dispersion curves are associated with specific values of 5, as shown in the panel legend. Generally, 
a finite gap separates the energies of the two bands; only at the two points 5 = = n /2 and <5 = ^2 = 37r/2 (blue and 

purple dashed lines, respectively) the gap vanishes for k — 0 and fc = tt, respectively, causing the presence of a discontinuity 
in the group velocity [see panel c)]. As shown in panel d), points 5i and 82 represent the boundaries between two topological 
phases, d) Topology of the Bloch eigenstates. States \(f>a{k)}, representing the polarization (coin) part of the eigenstates of the 
system, can be represented as points on the surface of a Poincare sphere, where they are individuated by ±ns{k). Here we 
represent these states as solid spheres, whose radius is related to the value of fc, as shown in the underlying legend. As the 
quasi-momentum k varies in the Brillouin zone {— tt, tt}, chiral symmetry forces these states to lie on a great circle of the sphere. 
The associated winding number W depends on 5, and determines the existence of a non-trivial and a trivial topological phases; 
the former {W = 1) occurs when < 8 < S 2 , while the latter {W = 0) is obtained in the remaining part of the interval {0, 27r}. 
As an example, from the top to the bottom of the panel, we represent the coin eigenstates of a QW with 8 — tt, 77r/4, 
respectively. For 8 = n these form a closed loop along the great circle, while in the other cases they go back and forth along a 
finite arch, whose lengths depends on the value of 8 (when 8 = {0, 27r}, the arch results in a single point). 
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FIG. 2. Layout of the setup. In the heralded single photon source (the optical components that are used to realize this 
source are not shown in the figure, for details see Ref. [19]) laser pulses (100 fs) at 800 nm generated by a Titanium-Sapphire 
source (Ti:Sa) with repetition rate 82 MHz shine a type I Beta-Barium Borate crystal (BBOl) for second harmonic generation 
(SHG); frequency doubled pulses at 400 nm, with 110 mW average power and linear-horizontal polarization, pump a type 
II BBO crystal (BB02), cut for collinear and degenerate spontaneous parametric down conversion (SPDG). Signal and idler 
photons, generated in horizontal and vertical linear polarizations, respectively, are spatially separated by means of a polarizing 
beam splitter (PBS) and then coupled to single-mode optical fibers (SMFs); the idler photon is directly sent to an avalanche 
photodiode (APD Dl, not shown) while the signal one, after passing trough the QW system, is analyzed in polarization and 
OAM and finally detected by APD D2, in coincidence with Dl. Before the QW, the photon exits the fiber in the 0AM state 
m = 0, and then its polarization is prepared in the state |<(>o) = q\L) -\- /3|i?); the two complex coefficients a and j3 (with 
|ap -I- |/3p = 1) are selected by using a half-wave plate (HWP) and a quarter-wave plate (QWP) (apart from an unimportant 
global phase). After the initial state preparation, the photon goes through the 6-step QW, with the single step consisting of a 
q-plate and a QWP oriented at 90°. For each g-plate, the value of the optical retardation 5 is controlled by the amplitude of 
an alternating electric field, which is introduced by means of an external generator [22] . At the exit of the QW, a polarization 
projection is realized using a second HWP-QWP set followed by a linear polarizer (LP). The OAM state is then analyzed by 
diffraction on a spatial light modulator (SLM), followed by coupling into a SMF. Before detection, interferential filters (IF) 
centered at 800 nm and with a bandwidth of 3.6 nm are used for spectral cleaning. The latter was required since the photons 
wavelength strongly influence the action of the devices implementing the QW (g-plate and QWP). 
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a) b) 




c) d) 




FIG. 3. Experimental results, a) Example of a probability distribution for the walker. Measured (green, left) and expected (red, 
right) probability distributions after a 6 steps QW of a photon initially prepared in the state m = 0 and {a, /3} = l/%/2 {1,1}, 
when S = 2.95. The error bars represent statistical errors at one standard deviation, calculated assuming Poissonian fluctuations 
on single counts, b-c) Measured values of statistical moments Mi/n and respectively, when varying S in the range 

{tt/SjTt} with steps of tt/IG, for initial states {a,/?} = 1/{1,1} (blue dots), and {q,/ 3} = {0,1} (purple dots). Dashed 
lines represent expected values for the plotted quantities, as obtained from numerical simulations. The continuous lines are 
the asymptotic limits reported in Eq. d) Measured values of y/M^jn when the initial polarization corresponds to the 
state cos(0/2)|L) + sin(0/2)|i?), varying the polar angle 9 G [0,7r] with steps of 7r/22; data are collected in correspondence of 
5 = 7r/4 (blue), 5 = Stt/S (purple), 5 — Srr/d (yellow), 5 — n (green). For each of these configurations, continuous lines give the 
corresponding asymptotic limit obtained using Eq.|^ When <5 = Stt/S, experimental data are slightly higher then the associated 
continuous line; nevertheless, the discrepancy is compatible with the results of numerical simulations for a hnite number of 
steps (see Fig. W)- These data show that M 2 is asymptotically independent of the input coin-polarization state. 

































SUPPLEMENTARY INFORMATION 


Appendix A: Role of the external parameter S in tailoring the structure of the energy band of the QW 

system. 


In this section, we discuss the dependance of the dispersion relations of our QW system with respect to the value 
of the parameter S. The explicit expression of the QW quasi-energies is given by the following relation 


Es{k) = ±cos 


cos {S/2) + sin {S/2) cos {k) 

71 


(Al) 


We recall that Es{k) are the eigenvalues of the effective Hamiltonian HeS', accordingly, they represent a quasi-energy 
defined in a Brillouin zone {—7r,7r}, as a consequence of the temporal coordinate (the step-number) being a discrete 
variable. The expression of the related group velocity, that is Vs = dEs/dk, is given by 


Vs{k) 


dEs{k) sin((5/2) cos(/c) 

dk a/2 — [cos((5/2) -I- sin(5/2) sin(A:)]^ 


(A2) 


In Fig.[^ we plot the quasi-energy and the group velocity as a function of the quasi-momentum k, for different values 
of <5 in both the trivial and the non-trivial phases. Here we restrict our attention to the upper band (the energy 
of the other has only opposite sign), and to half of the Brillouin zone, having the dispersion relation the symmetry 
E{k) = E{—k). At first glance, it is clear that the energy bands have marked differences with respect to the system 
being in the trivial or the non-trivial topological phases. In panel a), we can observe that as S goes from 0 to 7r/2 
(or from 37r/2 to 27r), the band covers a larger (smaller) range of energy values; accordingly, in panel c), we can 
observe that the maximum values for the group velocity increases (decreases). In contrast with what we observed 
in the trivial phase, in panel b) it can be noted that in the non-trivial phase a change of S results in a shift of the 
whole band, plus a tiny deformation of the curve. In panel d), we can observe that this results in a alteration of the 
group velocity dispersion that keeps constant the associated maximum value; this can be seen as a consequence of the 
geometric properties of the unit vector ns{k) in the non-trivial phase, being valid the relation Vs = = —riy. The 

consequences of these features on the statistical moments associated with a QW evolution in the large step-number 
limit have been deeply discussed in the main text. 


Appendix B: SSH model and QW system: analogies and differences 

A typical example of a ID periodic system with an internal 2-state degree of freedom, characterized by chiral 
symmetry, is the Su-Schrieffer-Heeger model (SSH) describing the electron dynamics in the polyacetylene chain. We 
discuss explicitly this model to introduce the symmetries that characterize our QW, highlighting the analogies and the 
differences between these different systems. The SSH model describes the polyacetylene chain as a one-dimensional 
(ID) dimerized lattice [7], where each unit cell consists of two sites A and B, as shown in Fig. la. Complex amplitudes 
t and t' quantify the hopping between adjacent sites, belonging to the same unit cell or not, respectively (for the 
purposes of our analysis it is sufficient to consider t and t' as real parameters). The Hamiltonian describing the electron 
dynamics along this dimerized lattice is typically introduced for a multi-particle system; accordingly we express it 
using the second quantization formalism, (please note that we adopted the typical choice to present the QW as a 
single particle model; in this context the two representations are completely equivalent): 

{ N N-1 'j 

^ ^ ^li,n+VB,n 1 ' -fh.C.. (Bl) 

n=l n=l J 

Here, „ and are the creation and annihilation operators for an electron in the cell n, on the sublattice A 

{B), and h.c. stands for “hermitian conjugate”. We are considering a finite chain made of N cells, which are spaced 
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a) c) 




b) d) 



FIG. 4. The parameter 5 and the dispersion relations of our QW system. We plot the dispersion relation for the upper energy 
band, and the associated dispersion of the group velocity, of a QW system described by the single step operator Uo, with 
varying the external parameter S in the range {0, 27r} with steps of tt/IO; in these plots, we restrict our attention to half of 
the Brillouin zone, being the dispersion relation symmetric in the other region {E{k) = E{—k)). In all images, coloured blue 
and purple curves correspond to the two configurations in which the quantum transition occurs, that is 5 = tt/2 and S = 37r/2, 
respectively. In panels a-b), coloured orange and green curves correspond to S being in the intervals {0,7r/2} and {37r/2,27r}, 
respectively. As a consequence of the energy bands properties, in the former case, as S increases, the maximum of the group 
velocity (panel b) becomes larger; in the latter configuration, this behavior is reversed. In panel c-d), coloured red curves 
correspond to <5 £ {7r/2, 37r/2}. In this case, the group velocity maximum remains locked to a constant value (panel d), while 
the corresponding quasi-momentum is drifting. In all panels, the correspondence between values of S and the associated plots 
reflects the order in which these curves are displayed. 


by a distance a. Using the Fourier transform for the creation and annihilation operators 

= (B2) 

^ n 

= (B3) 

^ n 

we can express the Hamiltonian operator in momentum space: 

H = ^ ipliit + t' cos (a k)] ax + t' sin (a k) ayy^jik- (B4) 

k 

In the latter equation, tpk is a 2D vector operator defined as 


= 



(B5) 


and ax and ay are the 2D Pauli operators. For a better comparison with the QW system, it is convenient to operate 
a shift of the quasi-momentum in the Brillouin zone, that is fc —> A: -|- tt. After this transformation, the eigenvalues of 
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the SSH Hamiltonian (B4) are given by: 


Et,t'{k) = — 2tt' cos (a k)]. 


(B6) 


Henceforth, we will assume a = 1, so as to switch to adimensional units (that characterize intrinsically the quantum 
walk system). Using Eq. B6 we can determine the group velocity dispersion, which has the following expression; 


= ±' 


tt' sin k 


\J — 2tt' cos (a k)] 

In momentum representation, the Hamiltonian can be conveniently expressed as a 2 x 2 matrix: 

H{k) = Et^t' (k) Tit^t' (fc) • cr 


(B7) 


(B8) 


where cr is the 3D vector whose components correspond to the three Pauli matrices. At each momentum k, the 
expression of the real 3D unit vector that determines the position of the coin eigenstates on the Poincar sphere 
representing the spin Hilbert space is: 


nt,t'{k) = 


t — t' cos k 


t' sin k 


t ^ — 2tt' cos k ’ — 2tt' cos k 


,0 


(B9) 


In the latter expression, we can observe that the ^ component of nt^t'{k) is vanishing; this results from the absence 
oi A — A and B — B links in the Hamiltonian. Accordingly, the vector nt^t'(k) is confined on the equator of the 
Poincare sphere that we use to represent the 2D sublattice Hilbert space. Moreover, we can note that the components 
of nt^t'{k) are related to the group velocity dispersion; indeed we have that 




riy = 


(BIO) 


I n Fig . we report the dispersion relations and the coin eigenstates for the SSH model, computed using Eqs. B6 
and |B9[ respectively, for different values of t and tk In these plots we can see that the two bands have a finite gap, 
vanishing when t = t' (Fig. Configurations t < t' and t > t' are not equivalent, being the winding number of 
rit^t'ik), as k varies in the Brillouin zone {—7r,7r}, equal to one (Fig. [^) and zero (Fig. [^), respectively; this is the 
same topological invariant that we introduced for the QW. Similarly to the latter system, we can observe that the 
dispersion of the group velocity is different in the two distinct phases; the main features are the same presented in 
the previous section, in the context of a QW evolution. 

Now, let us discuss briefly the symmetries characterizing the SSH system. First of all, the model has time-reversal 
symmetry: THT~^ = H, where the time-reversal operator T is defined by T = AT and K is the complex conjugation 
operator. Indeed it is straightforward to show that H*{k) = H{—k). Moreover the system has the sublattice, or chiral, 
symmetry SH{k)S~^ = —H{k). The sublattice symmetry operator S is defined by where is the z-component 
of the vector that represents the Pauli matrices and acts on the sub lattice degree of freedom. Finally the system 
has the particle-hole symmetry: PH{k)P = —H{k), where the particle-hole operator is defined by P = a^K. Then 
the topological class is BDI since P^ = = 1 [2^26) . The class of ID Hamiltonian operators with these 

symmetries is generally labeled as SSH, since the associated physical systems have the same topological features as 
those emerging in the SSH model. This is the case of the QW system we introduced, as we will demonstrate in the 
following paragraphs. 

As explained in the main text, a QW consists in the dynamical evolution of a system made of a walker and a coin part; 
in our photonic platform, they are encoded in the spin angular momentum and in the orbital angular momentum of 
light, respectively. In particular, the unitary operator Uq defining one step of the quantum walk, consists of a quarter 
wave plate, oriented at 90° with respect to the horizontal direction, followed by a g-plate. 


Uo = QsW, 


qwp • 


(Bll) 


The QWP action is described by the operator lUqwp, acting only on the polarization degree of freedom: 

IUqwp|T,m.) = ^ {\L,m) + i \R,m)) 

—t 

Wqwp|P, m) = ^ (|L, m) - i |P, m)) 


(B12) 












11 



FIG. 5. Dispersion relations and eigenstates for the SSH model, a-b) The values of the quasi-energy for the two bands are 
reported as a function of the quasi-momentum k, in the Brillouin zone {—rr, tt}. Coloured red and green curves correspond to 
the two configurations = (1,1.5) and = (1,0.3), respectively. The dashed curve corresponds to the case t = t' = 1. 
b) Group velocity dispersion for the lower energy band. Colour legend is the same with respect to panel a). From these plots, 
we can observe that the energy gap closure at k — 0 (dashed curves) is associated with a jump of the group velocity, c-d) 
Representation of the vector nt^t'(k), calculated when varying k in the range {—7r,7r}, with steps of 7r/4; hopping amplitudes 
(t, t') are equal to (2,1) and (1, 2), respectively. These vectors are represented by solid spheres, whose radius is associated with 
the quasi-momentum value, as reported in the legend. The coin eigenstates corresponding to the two configurations t < t' 
[panel c)] and t > t' [panel d)] have marked differences in terms of their topology. Indeed, when varying k in the Brillouin 
zone, the associated trajectory over the Poincare sphere results in closed loop or in a finite arch, respectively. Accordingly, the 
winding number is one in the first case, whereas it vanishes in the second one. 


The g-plate action is described by the operator Qs 

Qs\L, m) = cos {5/2)\L, m) + i sin {S/2)\R, m + 2q), 

Qs\R,m) = cos (5/2)1 ii, m) -|-t sin (5/2)|L, m — 2q). (B13) 


Here, q is the topological charge of the g-plate; \L) and |i?) represent left circular and right circular polarizations; 
and 6 denotes the optical birefringent phase-retardation. In our experiment, we considered q = 1/2. It is possible 
to construct a time independent effective Hamiltonian, Hes, such that the quantum walk represents a stroboscopic 
simulator of the evolution generated by Hgg at discrete times. By using momentum eigenstates jfc), iJeff is given by: 


^eff(5) = dk [Es{k)ns{k) ■ a] 0 \k){k\, 


(B14) 


where Es{k) represents the modulus of the Hamiltonian eigenvalues (see Eq. Al), cr denotes the vector of Pauli 
matrices, and ns{k) is a real 3D unit vector representing the coin eigenstates on the Poincare sphere associated with 
the spin Hilbert space. The components of this vector are given by 


nxik) = {cos (5/2) — sin (5/2) cosfc} /N{k), 
^yi^) = — sin(5/2)sinA:/7V(fc), 

nz{k) = -ny{k), 


(B15) 
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where the quantity TV is a normalization factor 


Nik) = ^2{l-cos^Esik)]}. 

It is worth to be noted that comparing Eq. |A2|and Eqs. |B15 it is possible to verify that 


Vsik) = = -r 


(B16) 


(B17) 


This important result, which links some geometrical aspects of the QW system (that are used to define the system 
topology) and the dispersion relation, has been discussed widely in the main text. It is straightforward to show that 
there is a vector a perpendicular to nsik), for all values of k and 5; in other words, these eigenstates are positioned 
on a great circle of the Poincare sphere, and such circle is independent of the value of S (see Fig. [^. For a simpler 
comparison between the QW system and SSH model, we can choose a as the quantization axis for the coin {z 
component). In the new representation, the single step operator Uq is obtained by means of a rotation of 7r/4 around 
the x-axis of the Poincare sphere in the coin space. Accordingly the coin part of the associated eigenstates is confined 
on the equator of such sphere, that is n(,(/c) = 0. It is straightforward to show that the effective Hamiltonian operator 
(B14), with ns{k) replaced by nj(fc), has the same symmetries introduced for the SSH model, thus it belongs to the 
same topological class. 


The SSH and QW systems we introduced have the same topological classification and may show distinct phases, 
according to the values of some external parameters. In SSH model, the quantum transition occurs when the hopping 
amplitude between different cells is equal to the one between the two sites of the same cell, that is t = t'. In our 
quantum walk system, the phase transition occurs at 5 = 7r/2 and 5 = iTTl2, when sin | = ±cos |, and the walker 
has the same probability either to remain in the same lattice site or to pass to the neighbour one. 

It is straightforward to show that our QW model, when the localized Wannier wavefunction basis is adopted, can 
be mapped into a SSH-like model defined on a one-dimensional lattice where, in addition to the link between the 
two sites within the same cell, only hoppings between nonequivalent sites of different cells are present. However, 
unlike the SSH model, in the effective Hamiltonian not only nearest neighbor cells connected: the hopping between 
nonadjacent cells is possible, that is the photon transfer integral associated with the orbital angular momentum 
displays a long-range behavior. This emerges when analyzing the property of the effective Hamiltonian, even though 
in such discrete-time systems the physically meaningful quantity to be considered is the single step evolution operator, 
which shows couplings only between nearest-neighbor cells. 


Appendix C: First and second order moments of the walker probability distribution: the large step-number 

limit for a QW system 


The large step-number limit for the walker distribution moments (Eqs. I][2 ) can be derived conveniently by evaluating 
Ml and M 2 in momentum representation, where they are defined as follows: 


Ml = 


dk 

dk 




M 2 = I ^ (</>o| (j7o^) 


(Cl) 


Here n is the step-number, Uq is the single step evolution operator (Bll), and |:/)o) is the coin initial state. Expanding 
the evolution operator as 


Uq = Exp{—inEsik) [nsik) ■ a]} 

= cos {n Esik)} I 2 — i sin {n Esik)} [nsik) ■ cr] 


(C2) 


where I 2 is the identity matrix in 2D, and the components of nsik) are those reported in Eq. B15 it is straightforward 
to obtain the following equations: 


dA- 

— Vsik) (())o| insik) ■ a) [(jjo) + 0(l/n), 
M 2 /n^= f ^ [Vsik)f + Oil/n^), 


(C3) 

(C4) 

(C5) 
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FIG. 6. Analysis of statistical moments |Mi, M 2 } in a QW. The system, initially prepared in the state |'i/'o) = |0o)c <8) |0)u,, 
undergoes a QW described by the step operator (Bll 1. For every plot, purple points are obtained from a numerical simulation, 
when varying S with steps of 7r/16 in the range {0, continuous blue lines represent the quantity L{S) [panels a-c)], or Ij{5) 
[panels d-f)] (see Eq. C91. For the simulation, we prepared the coin in the state {si, S2, S3} = {0,0,1}. Panels a-c) First order 
moment Mi, divided by the number of steps of the walk, as a function of the parameter S, for a walk of 5, 20, and 50 steps, 
respectively (this number is reported inside each figure). As n increases, simulated data converge to the values predicted by 
Eq. 0 Panels d-f) Square root of the second order moment, divided by the number of steps n. The figures are organized as in 
panels a-c). In this case, we can observe that simulated data converge much faster to the asymptotic values reported in Eq. 
with a discontinuity emerging even for a walk of few steps. 


Eq. 


BIT) 


C3 reporting the expression for the first order moment, can be evaluated considering the Vs = Uz = —ny (Eq. 
Applying this substitution, the same equations reads 


dA- 

^ [ys{k)f + 0{lln) 


(C6) 


A plot for L[5) (and for sfL) is reported in Fig. ^ where this asymptotic limit is compared to simulated data for a 
n-step QW, with n = 5, 20, 50. As discussed in the main text, it can be observed that as n increases, both Mi and 
M 2 converge to their asymptotical limit with this convergence being much faster for M 2 with respect to Mi. 

First and second order moments reported in Eqs. |C6|C4| coincide with the expressions reported in the main text 
(Eqs. |l|2[ ). In the specihc case of our QW model, the integral appearing in both Eqs. C6|C4 can be solved analytically. 
Using the expression for the group velocity reported in Eq. |A2[ Eq. [^yields 




^ 271 - 


m = 


2tt 


Vidk = 


sin^((5/2) cos^(fc) 


2tt Jq 2 — [cos((5/2) -I- sin(5/2) sin(A:)]' 


-dk. 


(C7) 


This integral can be calculated by the residue theorem passing to the complex variable z = e*^. Then we obtain 
L{S) = j> fs{z)dz, where the integral is along the unit circle in the complex x-plane and fsiz:) is given by 


fsiz) = 


i(l -I- z^)^ sin ((5/2) 


7rz[(l -f cos(5) — — lOz^ — 1 — 4fz(z^ — 1) sin((5)] 


(C8) 


The poles of fs{z) are located on the imaginary axis at Zk = — 1) cot(i5/4), f(-\/2 -|- 1) tan(i5/4), —i(-\/2 -|- 

1) cot(5/4), —i{y/2 — l) tan(5/4)} (fc = 1,..., 5) and the residues of fs{z) at the poles are given by 27rirfc = {1, \\—\/2 + 
2cos(5/2)], i[v^- 2cos((5/2)], i[y2-b2cos((5/2)], \[-^/2- 2cos((5/2)]}, respectively. Apart from the pole at z = 0, 
when 5 varies, the locations of the poles move in the complex plane entering and exiting the unit circle, but only the 
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residues of the poles inside the unit circle contribute to L{S). Thus, we find 


L{6) 


27ri(ri + ra + ra) = 2sin^((5/4) 

' 27ri(ri + r2 + rs) = 1 - 

27ri(ri + r 2 + r 4 ) = 2 cos^(5/4) 


for 0 < (5 < 7 r /2 
for 7 r /2 < (5 < 37 r /2 
for 37 r /2 < (5 < 27r 


(C9) 


Appendix D: Analysis of first and second order moments for the SSH model 

The same approach we used for the analysis of statistical moments in the QW system can be used for the SSH 
model. In this case, we have that 


Mijr = -S 2 £{t, t') -b 0 (l/r), 

(Dl) 

M2/t^ = £(£ t') + 0 ( 1 /t^), 

(D2) 


where we are considering t as a continuous temporal coordinate. The quantity C has the same expression reported 
in Eq. with the group velocity being that of the SSH model (Eq. B7). It can be noted that these expressions are 
obtained when considering an electron prepared in a localized initial state. While this is a standard choice for QWs, 
in electron dynamics this condition is hard to be reproduced experimentally, and typically is not considered. Even in 
this case, C has an analytical expression. Passing to the complex variable z = e®^, we find 


C = - 


i tt' 


dz 


(z^ - 1)= 


StT / 2;2 [^2 - Z + l] 


(D3) 


where the integral is evaluated along the unit circle |z| = 1 in the complex plane. Eq. D3 can be solved by residue 
method. The integrand function has three poles, which are located along the real axis: 


zo = 0 ; zi = t/t'; Z 2 = t'/t, 


(D4) 


It is worth noticing that zq is doubly degenerate. When varying t and t', poles zi and Z 2 move in the complex plane, 
entering or exiting the unit circle. Since only inner poles contribute to the integral ( |D3[ ), residues at zi and Z 2 should 
not be considered simultaneously; thus, in the expression of £ the residue at zi (Z 2 ) will appear when t < t' {t > t'). 
If t = P, zi and Z 2 are located on the integration path (the unit circle) and the Eq. |D3| cannot be solved using this 
method. The residues at the poles z : k reported in Eq. |D4| are 


27riro = 


r +t 
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27ri ri = 


£ - t 


,/2 


2'Kir2 = 




(D5) 


Accordingly, we have that 


£{£t') = 


27ri(ro + r2) = ^ 
27ri(r-o +ri) = ^ 


for t' <t 
for t' > t 


(D 6 ) 


In Fig. [^we plot the function £, and we compare it to simulated results for the evolution of the electron dynamics 
described by the SSH model. In the large time limit, numerical data converge rapidly to £ (D 6 I. It is important to 
observe that, as a difference with respect to the QW model, here M 2 would vary in the non-trivial phase too, for any 
variation of the parameter t, that we considered as a constant. The observable quantity that is locked to a constant in 
the non-trivial phase, independently of the hopping parametrization, is M 2 /t^. Remarkably, as for the QW system, 


this is equal to Uy, integrated over the Brillouin zone in momentum space (see Eq. BIO). 
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1.0 1.5 2.0 


t- (t=l) 


FIG. 7. Second order moment for the SSH model in the large time limit. We plot the second order moments M 2 corresponding 
to the probability distribution for a single electron whose dynamics is ruled by the SSH Hamiltonian (B11. As for the QW, we 
considered a localized initial state (at r = 0). Blue circles represent the data obtained from a numerical simulation, in which we 
evolved the initial wave function at the time r = 50 (we recall that we are dealing with adimensional units). In such simulation, 
we considered t as a constant parameter (i.e. t = 1), while varying t' within the interval {0,2} with steps of 2/25. The red 
continuous line represents the asymptotic limit calculated using Eqs. |D2| and |D6| We can see that such limit reproduces well 
the simulated data. 








